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Abstract 

A synchronization mechanism driven by annealed noise is studied for two rephcas 
of a coupled-map lattice which exhibits stable chaos (SC), i.e. irregular behavior 
despite a negative Lyapunov spectrum. We show that the observed synchronization 
transition, on changing the strength of the stochastic coupling between replicas, be- 
longs to the directed percolation universality class. This result is consistent with the 
behavior of chaotic deterministic cellular automata (DC A), supporting the equiv- 
alence Ansatz between SC models and DCA. The coupling threshold above which 
the two system replicas synchronize is strictly related to the propagation velocity 
of perturbations in the system. 
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1 Introduction 



The occurrence of disordered patterns and their propagation in the presence 
of a negative Lyapunov spectrum have been often observed in spatiotemporal 
systems [1-5] . One can classify this kind of irregular behavior into two general 
groups with basically different features: transient chaos and stable chaos. 

Transient chaos (TC) is a truly chaotic regime with finite hfetime, character- 
ized by the coexistence in the phase space of stable attractors and chaotic 
non attracting sets - named chaotic saddles or repellers [6]. The system, start- 
ing from a generic configuration, typically exhibits irregular behavior until it 
collapses abruptly onto a nonchaotic attr actor. 

Stable chaos (SC) constitutes a different kind of transient irregular behavior 
[1,2] which cannot be ascribed to the presence of chaotic saddles and there- 
fore to divergence of nearby trajectories. In SC systems, moreover, the time 
spent in transient regimes may scale exponentially with the system size (su- 
pertransients [1,2]), and the final stable attractor is practically never reached 
for large enough systems. One is thus allowed to assume that such transients 
may be of substantial experimental interest and become the only physically 
relevant states in the thermodynamic limit. While TC remains associated to 
information production, i.e. to the response of the system to infinitesimal dis- 
turbances, SC is mainly related to propagation and mixing of information. In 
other words, SC systems are sensitive only to perturbations of finite ampli- 
tude [7], while they respond to infinitesimal disturbances in a way similar to 
stable systems. Such a key feature makes meaningless any characterization of 
the SC complexity by means of Lyapunov theory. 

In this paper we focus on SC behavior whose origin has not yet found a con- 
vincing explanation, even though, it has been observed in several spatially- 
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extended models such as coupled map lattices [1-3] and oscillators [4]. To 
provide an adequate description of SC systems we invoke their strict simi- 
larity with discrete models such as deterministic cellular automata (DCA). 
In fact, according to the conjecture that "SC systems represent a continuous 
generalization of deterministic cellular automata" [2], we can argue that, what 
is known about DCA could be automatically translated into the language of 
SC. Although a general mapping of SC- onto DCA-models is still missing, 
the conjecture is supported by the fact that, also in finite size DCA, limit 
cycles and fixed points are the only allowed attractors, since the number of 
possible configurations is finite. Moreover in some DCA, the transient dynam- 
ics towards the final attractor may exhibits a long living irregular behavior 
with lifetimes that typically grow exponentially with the system size. In this 
case it is practically impossible to find any recurrence (Poincare cycles) for 
large systems. According to Wolfram classification [8], these DCA form the 
third ("chaotic") class and they share several properties with continuous SC 
systems. 

The emergence of this "chaoticity" in DCA dynamics, is effectively illustrated 
by the damage spreading analysis [9,10], which measures the sensitivity to 
initial conditions and for this reason is considered as the natural extension 
of the Lyapunov technique to discrete systems. In this method, indeed, one 
monitors the behavior of the distance between two replicas of the system 
evolving from slightly different initial conditions. The dynamics is considered 
unstable and the DCA is said chaotic, whenever a small initial difference 
between replicas spreads through the whole system. On the contrary, if the 
initial difference eventually freezes or disappears, the DCA is considered non 
chaotic. 

A similar scenario holds for systems exhibiting SC [3], where a transitions 
between regular and irregular dynamics may occur upon changing a control 
parameter (e.g. the spatial coupling between sites); the irregular dynamics is 



3 



often associated to spreading of damages. 



In this work we show that another characterization of the SC behavior can be 
achieved through a suitable synchronization method, which has been success- 
fully employed in Ref. [11] to classify the chaotic properties of DCA. In our 
opinion this method may be used to complement the common damage spread- 
ing analysis. The basic idea consists of measuring the minimal "strength" of 
the coupling between a replica (slave) and the original system (master) re- 
quired to achieve their perfect synchronization. The master-slave interaction 
is a stochastic and spatially-extended version of the Pecora- Carroll synchro- 
nization mechanism [12]. Unlike the damage spreading method, where replicas 
are independent, this coupling scheme implies that the evolution of slave is 
driven by the master. Practically each time step is composed by two phases: in 
the first one, the master and slave system evolve freely with the same evolution 
equation, and then a fraction p of the degrees of freedom in the slave-system 
is enforced to take the value of the corresponding degrees of freedom in the 
master. Synchronization of spatially extended systems has been usually stud- 
ied with symmetrical interactions [13-17]. Our asymmetric scheme, instead, 
allows probing the dynamical properties of the master system. Through a 
gradual increase of p from to 1, the dynamics of the slave system tends to 
synchronize to that of master, and at a threshold p* a synchronization tran- 
sition occurs: the pinching synchronization transition (PST). The threshold 
p* above which the replicas synchronize, is an indicator of the chaotic behav- 
ior of the unperturbed system (master). Indeed, a large value of p* implies a 
large fraction of sites to be pinched in order to achieve the synchronization, 
indicating that the dynamics of the replicas is rather irregular and difficult to 
control. 

In DCA, the PST belongs to the directed percolation (DP) universality class [11,17], 
while in fully chaotic continuous systems (e.g. coupled map lattices) the syn- 
chronization is never perfect for finite times, nor it is equivalent to an absorbing 
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state [17]. This generally implies non-DP scaling exponents [18-20]. Here, we 
find that the PST is well defined for a model showing typical SC behavior 
originally introduced in Ref. [2]. This PST is found to belong to the DP uni- 
versality class, in agreement to what happens for DCAs. The result further 
supports the conjecture that SC class contains DCAs. 

The sketch of this paper is the following. In the next section, we study a simple 
coupled-map lattice which is particularly suitable for discussing the distinction 
between SC and TC behaviour. The model, in fact, displays a TC regime before 
falling into DCA dynamics with the typical properties of a discrete SC regime. 
In Section III, we describe the SC model of Ref. [2], whose dynamics never 
reduces to a DCA, and the pinching synchronization technique applied to it. 
For the latter model we obtain the PST phase diagram (Sec. IV), by measuring 
the synchronization threshold as a function of the coupling strength between 
sites. Such a phase diagram agrees remarkably with that ah^eady obtained 
through damage spreading analysis in Ref. [3]. We provide an argument to 
explain this consistency. Finally, conclusions and remarks are reported in the 
last section. 



2 Transient and stable chaos 

It is instructive to discuss the main differences between SC and TC regimes 
with the aid of a simple spatiotemporal model in which both of them occur. 

Let us consider the one-dimensional coupled-map lattice (CML), i.e. an array 
of state variables {xi, ...,xl} in the interval [0, 1] subject to the discrete-time 
evolution rule 

Xi{t + 1) = ge{Xi^l{t),Xi{t),Xi+i{t)). (1) 
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Each site-variable Xi interacts diffusively with its nearest neighbors 



g,{u, V, z)^{l- 2e)f{v) + e [f{u) + f{z 
where e sets the coupling strength among maps. 

The local mapping has the form shown in Fig. 1: 

for < X < a, 

for q; < X < 1/3 — a, 



(2) 



X — a 



1 - 



1/3 -2q; 



X — 1/3 — a 
1/3 -2a 



X — 2/3 — a 
1/3 -2a 



for 1/3 — a < X < 1/3 + a, 

for 1/3 + a < a; < 2/3 -a, 

for 2/3-a<x <2/3 + a, 
for 2/3 + a < a; < 1 

for 1 — q; < X < 1. 



(3) 



a. 



In this example, we always use the democratic coupling £ = 1/3, and we shall 
neglect to indicate the e dependence. A typical grayscale pattern generated 
by the evolution of the above CML is shown in Fig. 2 (see the caption for the 
grayscale-code), for a = 0.068. Note that the continuous dynamics (gray) is 
limited to restricted domains, while in the rest, the system has fallen into a 
Boolean dynamics (typical of DCA), which is however far from being trivial. 

Systems showing SC and TC regimes are asymptotically stable, therefore the 
standard Lyapunov spectrum is not able to detect their transiently irregular 
states. In TC regimes, however, the existence of chaotic saddles in the phase 
space is revealed by the finite-time (or effective) Lyapunov exponent 



7(t) = ^(log||w(t)||), 



(4) 



where w(i) = {wi{t), ...,WL{t)} indicates a generic infinitesimal perturbation 
(i.e., a tangent vector) at time t, which evolves following the linearization 
of Eqs. (1-3). The average in expression (4) is taken over the ensemble of 
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trajectories which have not yet left a certain neighborhood of the saddle at 
time t [21]. The indicator ^(t) is expected to fluctuate around a positive value 
during the transient [5] and switches to negative values after the transition 
to the stable attractor occurs. In SC, however, even the finite-time Lyapunov 
exponent does not provide much information, as it becomes negative already 
in the transients. This indicates that the source of the SC behaviour cannot 
rely on the instabihty associated to repelling sets. Although, the existence of 
repellers may not be excluded a priori, certainly their role is not observable 
in SC. 

We now see how the above considerations apply to our toy-model. First, wc 
discuss the case of infinite slope map, (i.e. a = 1/6, full line in Fig. 1) corre- 
sponding to a pure SC regime. In fact, after one time step, each configuration 
of the lattice reduces to a sequence of "0" and "1" (Boolean configuration). 
The system evolution remains, however, irregular since, when a — 1/6, the 
model is equivalent to a DCA which follows the rule 150 

g{u, v,z) = u + v + z — 2{uv + uz + vz) + Auvz 

with u,v and z Boolean variables. This dynamics is known to generate highly 
irregular patterns [8]. On the other hand, this irregular behaviour cannot be 
associated to either chaotic saddles or local fiuctuations of Lyapunov exponent, 
the latter being — oo due to the specific form of the map. Despite Lyapunov 
analysis ensures that this system is totally insensitive to infinitesimal per- 
turbations, finite perturbations of amplitude greater than 1/6 give rise to a 
"defect" which propagates through the lattice. Indeed, it can be easily shown 
that a defect also evolves with the chaotic rule 150, because the rule is additive 
modulo two. 

A shght tilting of the vertical edges of the map, (0 < a < 1/6, dotted and 
dashed lines in Fig. 1) introduces some expanding regions in the phase space. 
Accordingly, one obtains a typical TC behaviour, due to the competition be- 
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tween stable and unstable effects, which decays into the above mentioned SC 
regime (see Fig. 2). Figure 3 shows a typical time fluctuation of the local ex- 
pansion rate (or local multiplier) fx{t) = ||w(t) ||/||w(t — 1)||, for a = 0.052, w 
being a tangent vector as in Eq. (eq:lyap). The irregularity of the signal is a 
clear indication of the chaotic-like behavior of the system. The simulation is 
stopped at time when w{t) becomes exactly 0, i.e., when the system set- 
tles down into a Boolean conflguration (SC state) which is Lyapunov-stable 
by deflnition. Therefore, provides an estimate of the time spent in the TC 
regime. 

For all those trajectories which have not yet entered binary configurations at 
time t, 7(t) remains positive, as seen in Fig. 4, where the distribution of ^{t) is 
shown for a set of 2000 trajectories starting from arbitrary initial conditions. 

The lifetime of the TC regimes preceding the SC behavior depends on a. 
becomes shorter as a approaches 1/6, value at which it vanishes, because the 
system behaves as a true DCA after just one time step. In the limit ck — > 0, 
instead, the flat regions of the map disappear (the slope being equal to 3) 
and TC regime becomes persistent and degenerates into fully developed chaos 
(FDC). 

A more detailed analysis of the behavior of the model of Eqs. (1-3) upon 
changing the control parameter a will be presented elsewhere. Here such a 
qualitative discussion has the only aim to highlight the key differences between 
SC and TC dynamics, both present in this toy model. 

In the example discussed so far, the SC regime occurs only as a DCA behavior, 
since, as soon as the system falls into a Boolean conflguration, it evolves as 
a genuine "chaotic" DCA. In the following section we discuss a continuous 
SC model whose behavior never reduces to DCA dynamics. In particular we 
study the synchronization properties of two replicas of such a system. 
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3 The model and the synchronization dynamics 



The dynamical system considered now, is the one-dimensional CML of Eqs. (1,2) 
with the coupling constant e e [0, 1/2] and periodic boundary conditions over 
a length L (system size) . The local mapping has the form 



as shown in Fig. 5. 

We use here the parameter values (a — 0.07, h — 2.70, c = 0.10) of Ref. [3] for 
which the map of Eq. (5) is attracted into a stable period-3 orbit. An example 
of the space-time evolution of the CML is shown in Fig. 6, where the presence 
of propagating structures similar to those in Fig. 2 is observed, despite here 
the system dynamics never relaxes onto a pure DCA state. 

In Ref. [3] an £-dependent dynamical phase transition between periodic and 
chaotic regimes of this system has been carefully investigated by damage 

spreading analysis. The periodic (chaotic) phase is characterized by the ab- 
sence (presence) of damage propagation which is found to behave linearly 



being S the linear size of the region affected by the damage. The factor 2 is a 
consequence of the symmetric coupling scheme which requires the left and right 
damage- front to progress at the same velocity Vp, but in opposite directions. 
The damage spreading velocity Vp can be considered a good indicator for this 
transition, because it vanishes in periodic phases. It turns out from Ref. [3], 
that for e < ~ 0.3 only periodic phases are observed, for e > ef^ ~ 0.3005 
only chaotic states exist. In the intermediate region e e ^i^^], periodic and 
chaotic behaviors alternate in an apparently irregular manner (fuzzy region). 



= < 



hx if < X < 1/6 

a + c{x-l/h) ifl/6<a;<l 



(5) 



Sit) = 5(0) + 2VFt 



(6) 
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We study the effects of the pinching synchronization on this model, and its 
interplay with the above described transition. 



The master system follows the dynamics of Eq. (1), while the slave system 
evolves as 

yi{t + 1) = [1 - ri{t)]g^{yi_i{t),yi{t),yi+i{t)) + ri{t)ge{xi_i{t),Xi{t),Xi+i{t)) 

(7) 



where ri{t) is a Boolean random variable 

Tiit) = < 



1 with probability p, 
otherwise. 



Practically, at each time step, a fraction p of site variables in the slave system 
is set equal to the corresponding variables of the master (see. Fig. 7). In the 
limit case p — 0, the slave system evolves independently of the master, while 
for p = 1, its evolution coincides with the master one. 

The synchronization order-parameter is the asymptotic value of the topolog- 
ical distance p between master and slave systems (i.e. the fraction of non 
synchronized sites) 

p{t,p)^]im Lj2^{\x,{t)-yM) (8) 

where ©(s) is the unitary step-function. We denote by p* the synchronization 
threshold, such that p(oo,p < p*) > and p(oo,p < p*) — 0. This synchro- 
nization mechanism defines an associated directed site-percolation problem in 
d= 1 + 1 dimension, where a site of coordinate (z,t) is "wet" if rj(t) = and 
it is connected to at least one neighboring wet site at time t — 1. At t = all 
sites are assumed to be wet. We denote by Pc the critical threshold for which 
a cluster of wet sites percolates along the time direction. The master and the 
slave systems can stay different only on the cluster of wet sites. 



10 



For chaotic CMLs two synchronization scenarios are possible, called weak and 
strong chaos respectively [22] . A system is said strongly chaotic if it does not 
synchronize even on the critical wet cluster (i.e. for p* = Pc) and therefore the 
distance p of Eq. (8) always exhibits DP scaling. Alternatively one can say that 
for strongly chaotic systems, the active and the wet clusters are essentially the 
same for every value of p. 

Instead, for weakly chaotic systems, the synchronization threshold p* is always 
located below Pc and the transition is discontinuous (first-order like) . For p < 
p* again the active and wet clusters coincide, whereas for p > p* the active 
cluster disappears, but the wet cluster still percolates. Such a behavior is due 
to the exponential vanishing of the difference field h{t) — {xiit) — yi{t)}f^i, 
even though local fluctuations of h can sporadically appear. 

However, for chaotic systems, the synchronized state is not robust with re- 
spect to an infinitesimal perturbation in the absence of the synchronization 
mechanism, i.e. it is not a proper absorbing state. Conversely, DCAs [11] al- 
ways synchronize at p* below p^ (i.e. synchronization occurs in presence of the 
percolating wet cluster) but p{t,p) still exhibits DP scaling. This is a straight- 
forward consequence of the finiteness of the number of states in DCAs, which 
prevents fiuctuations of h in the absorbing state. The toy model of Section II 
trivially follows this behavior, and the synchronized state is insensitive with 
respect to sufficiently small perturbations. We show in the following that this 
scenario holds also for continuous SC systems. 

4 Numerical results 

A first set of simulations has been performed for lattices of size L — 3000, 
with periodic boundary conditions. We measured p{t,p) at different times and 
the results have been averaged over a 5000 randomly chosen initial conditions. 
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For each simulation a transient of 10^ time steps has been discarded in order 
to avoid initial bias and reach stationary states. A site i is considered to be 
synchronized (and yi{t + l) is set equal to Xi{t+ 1)) if \yi{t) —Xi{t)\ < r, where 
r is a sensibility threshold. In this way we can control the effects of the finite 
precision of computer numbers. We checked that the results are independent 
of T (for small r) , and we choose r = 10~^ for massive simulations. 

A typical behavior of p{t,p) near p* is shown in Fig. 8. The value of p corre- 
sponding to the most straight curve at large t represents the best approxima- 
tion of p*, which can be estimated with good accuracy. The straight dashed 
line indicates the DP scaling results. 

For those values of £ for which the above analysis provided a too uncertain 
result in the estimation of p*, we have carried out single-site simulations for 
larger systems, obtaining a more accurate determination of p*. These simula- 
tions consist in preparing the slave system exactly synchronized to the master 
except for the central site, and in measuring how non-synchronized sites prop- 
agate throughout the system, generating DP-like clusters. This method probes 
the stability of synchronized states with respect to minimal perturbations. 

As usual in this type of simulations, we measured the survival probability P{t) 
of the desynchronized states, the number of non synchronized sites N(t), and 
their second moment R'^ (t) with respect to the center of the lattice (generally 
called the gyration radius) . Near the synchronization threshold and in the long 
time hmit, these magnitudes are expected to scale as [23,24] 



The determination of the asymptotic value of such exponents (e.g. 77) is per- 
formed by plotting the effective exponent 



N{t) - t", P{t) - t 



-s 



(9) 



log(iV(at)) 
log(iV(t)) 
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versus 1/t for several values of p. Here a indicates an arbitrary scale factor 
and we always set a = 2. In the limit i — > oo 77 converges to 77 for p — p* ■, 
and diverges for other values of p. In pure DP systems this method allows 
the simulation of effectively infinite lattices, since the reference state, i.e. the 
absorbing one (usually made of "0"s) is unique and does not change in time. 
Conversely, in our case the absorbing state coincides with the synchronized 
state and this imphes the detailed knowledge of the evolution of the mas- 
ter system. This circumstance imposes severe limitations on lattice sizes and 
performances of the method. 

In Figure 9 we report the behavior of the effective exponents rj, 5 and + 7] 
versus 1/t for several values of p. The asymptotic values of the exponents 
(77 = 0.330(5), 5 — 0.13(2), z — 1.25(2)) are consistent with the best known 
DP ones (77 = 0.31368(4), 5 = 0.15947(3), z = 1.26523(3) [24]). 

This analysis, repeated for several values of £, indicates that the synchro- 
nization of SC systems reasonably belongs to the DP universality class. The 
very slow convergence of the system dynamics to the asymptotic state makes 
hard to exclude that for some values of e the DP character of the transition 
is violated. However, we beheve that, owing to the stabihty of the model of 
Eqs. (1,2,3), small local disturbances are re-absorbed at exponential rate, and 
cannot generate desynchronization effects. In other words, the synchronized 
state is absorbing with respect to small fluctuations. This is also consistent 
with the fact that the simulation results are independent of the precision 
threshold r. Therefore, the local stabihty guarantees that, for what concerns 
synchronization properties, SC systems behave mainly like discrete ones and 
differently from continuous chaotic systems. This further supports the ^^equiv- 
alence Ansat£^ between SC and DC A. Indeed, preliminary simulations show 
that the DP scenario holds even if the local pinching is performed only up to a 
small difference A. This remark suggests the possibility of defining a finite-size 
maximum Lyapunov exponent, in a way similar to Ref. [25]. Further work in 
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this direction is in progress. 



Finally, the PST phase diagram is shown in Fig. 10, where p* is plotted versus 
e. The behavior of p* is compared with that of the damage spreading velocity 
Vp (Ref [3]), properly re-scaled. The consistency of the two phase diagrams, 
even in the fuzzy region, suggests that the indicators, Vp and p*, characterize 
different aspects of the same phenomenon. The strict correlation between Vp 
and p* is not surprising if one considers the changing rate of the density 
of non synchronized sites, n{t) — N{t)/L, in the pinching synchronization 
mechanism. For large time and systems and in a mean-field description of the 
process, wc can write for n{t) an equation similar to that one employed in 
contact processes and DP theory [26], 



The first term of the r.h.s. represents the active-site production due to the 
linear spreading of desynchronization regions, which occurs with a velocity Vp 
(see Eq. 6)). The second term accounts for the destruction of active sites by 
the pinching mechanism. The process becomes critical when the linear contri- 
butions, 2Vp n and pn balance, thus we obtain p* ~ Vp. The proportionality 
holds even outside of mean- field approach, i.e., when other terms, such as, 
higher powers of n, diffusion and multiplicative noise, are included in Eq. (10), 
provided one considers the renormalized parameters. The conclusion that p* 
is a nondecreasing function of Vp is however intuitive, because the higher Vp, 
the grater is the desynchronization rate, so the pinching probability has to be 
large to ensure the synchronization. 

5 Conclusions 

In summary, we have applied the pinching synchronization method to systems 
showing stable chaos. As for cellular automata, we have found that even in 




(10) 
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this continuous case the pinching synchronization transition (PST) is well 
defined, and that this transition belongs to the DP universality class. Our 
results show that the stable chaos is indeed equivalent to cellular automata 
"chaoticity" and definitively different from transient chaos. The PST phase 
diagram is consistent with that reported in Ref. [3] for damage spreading 
velocities, including the fuzzy region. 
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Figure captions 



Figure 1: Plot of the map of Eq. (3) for different values ol a: a — 0.02 (dotted 
line), a — 0.08 (dashed line) and a = 1/6 (continuous hne). 

Figure 2: Grayscale representation of space-time evolution of model of Eq. (2) with 
£ = 1/3 and a = 0.068. Time runs from top to bottom, white (black) 
color indicates Xiit) = (xj(t) = 1), while gray indicates all other values. 



Figure 3: Time behaviour of Lyapunov multiplier ^{t) — |z(i)|/|z(i — 1)| for CML 
of Eq. (2) with a — 0.052. The simulation is stopped at time T^. — 4455 
when the systems reaches the absorbing SC state. 

Figure 4: Histogram of the finite-time Lypunov exponent 7 of model of Eq. (2) 
computed from a set of 2000 arbitrary initial conditions, with a = 0.052. 
The size of the system is L = 3000 sites. For each initial condition, the 
evolution of Eq. (2) and its linearization are iterated untill the system 
approaches a Boolean configuration. 

Figure 5: Plot of the map of Eq. (5). 

Figure 6: Space-time evolution of the map Eq. 5, with e — 0.32. Grayscale from 
Xi{t) — (white) to Xi{t) — 1 (black). 

Figure 7: Schematic representation of the pinching synchronization. Vertical dashed 
lines indicate the sites of the slave system identified to those of the master 
corresponding to = 1 in Eq. (7). 

Figure 8: Log-Iogplot of the order parameter p{t,p) vs. t, forp = 0.1269,0.1270,0.1271. 
from top to down, and for e — 0.305. The dot-dashed straight line indi- 
cates the critical DP scaling t"^ with 5 — 0.159. The estimated synchro- 
nization threshold turns to be p* — 0.1272(1). 
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Figure 9: effective exponents ri{t) (a), S{t) (b) and z{t) + rj{t) (c) for e — 0.3004 
and several values of p. The average is taken over 50.000 runs, a — 2 and 
L = 2000. 

Figure 10: Dependence of the synchronization threshold p* on the coupling constant 
e (diamonds) compared with the behaviour of Vp (open circles), the Vp- 
values are properly re-scaled. In the inset an enlargement of the region 
around e — 0.3 is shown. The hue is a guide to eye. 
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